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■^ Anomalous superconformal Ward identities and commutator algebra in TV = 1 super- Yang-Mills 

P-i! theory give rise to constraints between the QCD special conformal anomalies of conformal com- 

^ posite operators. We evaluate the superconformal anomalies that appear in the product of renor- 

'-^ ' malized conformal operators and the trace anomaly in the supersymmetric spinor current and 

• ^ ! check the constraints at one-loop order. In this way we prove the universality of QCD conformal 

^ ■ 
^ anomalies, which define the non-diagonal part of the anomalous dimension matrix responsible for 
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■ ■ ■ ' scaling violations of exclusive QCD amplitudes at the next-to-leading order. 
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1 Introduction. 

Supersymmetric models [||], having higher space-time symmetry as compared to conventional 
ones, provide a strong consistency requirement on theoretical predictions. For the purposes of 
testing massless QCD calculations an especially illuminating example is A/" = 1 supersymmetric 
Yang-Mills theory, since both models have, up to a difference in colour representation of fermion 
fields, the same Lagrangian. Thus, we can map a QCD result to an A/" = 1 super- Yang-Mills 
theory one by identifying the colour Casimir operators in corresponding representations, i.e. Ca = 
Cp = 2NfTF = Nc- After this procedure a QCD result has to fulfill constraints arising from 
supersymmetry which connects gluonic and quark sectors of the theory. In this way the use of 
supersymmetry has allowed to find a set of identities p|, ^ between the entries of the forward 
anomalous dimension matrices of leading twist-two composite operators. They were valuable to 
clarify subtleties appearing in two-loop computations of anomalous dimensions. On tree level both 
theories are invariant under conformal transformations. Thus, the Af = 1 supersymmetric Yang- 
Mills theory is also invariant under superconformal transformations 0, |^ , which can give rise to 
a new set of constraints for certain conformal quantities which appear in the special conformal 
Ward identities for composite operators, the so-called special conformal matrix. Besides the 
breaking of conformal symmetry on quantum level by the trace anomaly in the energy-momentum 
tensor, we also have to deal with a superconformal anomaly due to the non- vanishing trace in a 
supersymmetric spinor current. Nevertheless, an explicit calculation of this anomaly will allow to 
check the special conformal anomalies calculated in QCD. 

Composite operators appear in various QCD applications by means of the operator product 
expansion and consequently their hadronic matrix elements contain a non-perturbative input, 
which is needed as an initial condition for the solution of the evolution equations. In the case 
of exclusive processes the off-forwardness of hadronic matrix elements, given in terms of distri- 
bution amplitudes and skewed parton distributions, requires operators with total derivatives. To 
ensure that the twist-two operators do not mix under renormalization at leading order in coupling 
constant, i.e. their anomalous dimension matrix has the diagonal form 'jjSjk, it is necessary to 
arrange the operators in such a way that they have a covariant behaviour under conformal trans- 
formations. This can be easily done. However, beyond one-loop approximation the anomalous 
dimension matrix jj^. develops non-zero non-diagonal, j > k, elements j^^ oc 0{al). 

The ordinary conformal algebra provides severe restrictions [^ 0] on the non-forward anomalous 
dimensions 7 of the conformal operators. In Refs. |]^, |^ we have developed a formalism based on 
the use of the broken conformal Ward identities for evaluation of the non-diagonal part, 7^°, of 
the complete anomalous dimensions matrix 7 = 7^ +7^°. This non-diagonal part arises entirely 
due to the violation of the special conformal symmetry in perturbation theory. The corresponding 



anomalies have been calculated to one-loop order accuracy in the minimal subtraction scheme 
using dimensional regularization, which imply the two-loop approximation for 7^°. To check our 
results, one can employ A/" = 1 super Yang-Mills constraints, valid in a renormalization scheme that 
respects supersymmetry, for the entries of the non-forward anomalous dimension matrix, derived 
in [Q. Unfortunately, this is not the case for the dimensionally regularized theory. Thus, one 
has to find finite renormalization constants from the latter to the dimensional reduction scheme, 
which is expected to preserve the supersymmetry. But there arises a subtlety in the evaluation 
of this rotation matrix for the gluonic sectoiQ which prevents it to be unambiguously fixed 0. 
Nevertheless, our result for two-loop non-forward anomalous dimensions is supported by the fact 
that the constraints can be fulfilled by a finite multiplicative renormalization, which proves the 
existence of a supersymmetric regularization scheme. 

Alternatively, we derive in this paper constraints directly for the special conformal anomalies at 
one-loop level and show that they are indeed satisfied. Our consequent presentation is organized as 
follows. In section |^ we define conformal operators, their anomalous dimensions, and relations of 
the latter to conformal anomalies. Section ^ is devoted to the study of translational and conformal 
super-anomalies on Lagrangian level in the dimensional regularization scheme. Then in section 
I we present transformation properties of conformal operators under relevant superconformal 
variations required for a derivation of the Ward identities discussed in the same section. In 
section |^ we give a derivation of relations between the scale and special conformal anomalies 
of conformal operators. Furthermore, we show that the latter acquire anomalous contributions 
originating from the product of the trace anomaly in the spinor current and conformal operators. 
They are explicitly evaluated in section |^, where it is demonstrated that indeed the anomalous 
constraints are fulfilled with special conformal anomalies from ^. Finally, we conclude. A few 
appendices are devoted to technical details that we found inappropriate to include in the body. 

2 Preliminaries. 

In this paper we discuss relations between the QCD scale and special conformal anomalies of 
conformal operators implied by the A/" = 1 supersymmetry. In A/" = 1 super Yang- Mills theory we 
introduce the conformal operators (for chiral even sector discussed throughout) 



«oi, = -r+itd^ycf ^ rr+, ^o], = G\%id^r'cf_\ |± t;^g^, (1) 



where CJ are the Gegenbauer polynomials and the tensor structures are r'^^'^) = (7+; 74.75), 
^llv'^^ = iajlv = 9tMu - n^nl - n*^ny] ie^upan^pria)- We use the convention d = d+d and V = V-V 
^This complication does not show up in the forward kinematics. 



with adjoint covariant derivative defined by Vl^ = dfj_S°''' + gf"''^''B'^. The '+'-sign as a subscript 
stands for contraction with the light-like vector n^ which specifies a direction along the light cone. 
For the latter purposes we introduce another vector n* such as n^ = n*"^ = and n ■ n* = 1. 
Obviously, the only difference from QCD arises in the gluino, which we loosely call quark, sector, 
which now belongs to the adjoint representation of the colour group. The factor | in Eq. (|ip is 
related to the Majorana nature of the quarks in the model. 

The renormalization group equation for these operators looks like 

T^[0,i] = -t^d^^l (2) 

where the square brackets will denote the renormalized operators defined by [Oji] = ^21=0 ^jkOki, 
with renormalization constant matrix Zjk, which generate finite Green functions with elementary 
field operators (p = {^; -S^}- We use everywhere the matrix notation and introduce the vector O = 
[gq] of quark and gluon conformal operators, which mix with each other under renormalization. 
As mentioned in the introduction, the only modification of a given QCD result is to identify the 
colour Casimir operators. Since at leading order the conformal anomalies have a unique colour 
structure^] it presents no difficulty to disentangle the separate components. 

We introduce as well the fermionic operator which is related to the bosonic ones (|I]) by super- 
symmetry 

^0;, = Gi'-^{zd^yp^''''>(^] r^r^. (3) 

Here P- ' are the Jacobi polynomials and the vertices read JF^^'"^) = (jnil'tl^)- The operators 
that form a representation of the supersymmetry algebra are defined by linear combinations of ([l|) 

= 0.;; «0;; + ^cuj ^O^^;, V,7 = ^/0^, U,i = g/0^, (4) 

with r = V{A) standing for the S (V) operator, and coefficients '^cu] = 1, '^to] = ^, ^to'^ = —^r, 
'^^'j = ^ &^d Qj = , j^t ■ Obviously, U = —75V. Note that the bosonic and fermionic 
conformal operators form the J\f = 1 chiral superfield 

$ = A + 2ex-e^J^, (5) 

with operators Sji and Vji {Sji and Vji) being the real and imaginary parts of the A (JF"!") complex 
fields, and Vj_i,/ identified with the Majorana fermion ( -. ) constructed from the Weyl spinor x- 

^More precisely, the anomalous dimensions in the gluon-gluon channel have in addition to the Ca term also 
trivial Nf dependent contributions, which arise from the self-energy insertion. 




Transformation between operators under supersymmetry arises from the conventional equation 
[CQ, ^]- = [CQ + QC ^]- = {Cr + fC} <?, with r = i^ and f = -2^ + 2^ ^. 

Now let us shortly point out, how the non-diagonal part of the anomalous dimension matrix 
is induced by the special conformal anomaly matrix. In four dimensional space-time the fifteen- 
parameter conformal group 5*0(4,2) is defined by its algebra containing the Poincare, dilatation 
T) and special conformal /C^ generators. The conformal anomalies are defined by the renormalized 
Ward identities. The generic form of the latter, however, written in an unrenormalized cast, reads 

{[Oji]5X) = -{5\0,i]X) - {t[0,i]5SX), (6) 

where X = JJ^ <P{xt) is a product of elementary fields appearing in the classical Lagrangian. 
Here 5 is any of the variations from the symmetry algebra in question. When the transforma- 
tion is a symmetry of the theory on quantum level then (JS* = up to possible BRST exact 
operators. In the (dimensionally) regularized theory the action does not vanish anymore for 
conformal, i.e. both scaling, (5'^0(x) = i[(l){x),T>\ = — {d^ + x^d^) (f){x), and special conformal, 
5^</)(a;) = i[4>{x), )Cfj] = — {2d^x^j_ — x^(9^ + 2xij_Xydij — 2iXuT,fj,y) (t>{x), variations, where d^ and S^,^ 
are the canonical dimension and the spin operator of the field 0, respectively. Thus, the renor- 
malization of the operator product i[Oji]6S is responsible for the conformal anomalies. 

Moreover, the commutator [P,/C_]_ = ilC^, where /C_ being the n*-light-cone projection of 
/C^, provides a connection between the conformal anomalies. In Ref. ^ the non-diagonal elements 
of the next-to-leading anomalous dimensions 

of QCD quark and gluon conformal operators were found in terms of one loop special conformal 
anomaly matrix 

aJ,HB)jf ^ -d,, (7?°^ - ^P^) + g,,. (8) 

It is constructed out of the leading order anomalous dimensions of conformal operators ^^^^^^ and 
the first expansion coefficient of the QCD /3-function (3q = |TpA^/ — yC^ times ajk{B)djf:{F)\ .^^ = 
—2{2k + 3) with ajk matrix from the conformal transformation of Ojk (see Eq. (^9]) below). 
The projector P = f ^ -^j in Eq. (H) singles out the gluonic component. Finally, the gr-matrix 
has appeared from the renormalization of the product of the conformal operator [Ojk] and the 
integrated trace anomaly S'^S oc — /(i'^x2x_6'^^(x) in the energy-momentum tensor 



j 



[0,l]5<lS = ^Y.Y,A0H-l] + ■■■. (9) 

fc=0 



In the dimensionally regularized theory, i.e. (i = 4 — 2e, the conformal variations of the QCD action 
can be calculated in a straightforward manner. Choosing the scaling dimensions of the physical 



fields equal to their canonical values in four dimensions^ {d^ = |, (i^ = 1) and setting the scaling 
dimensions of the ghost fields as d^j = and d^ = d — 2, the result 

5^8 = j d'^x w^(a;)| - ^ (Oa{x) + Ob{x) + Qa.{x) - Q^^{x) - Qd{x)) 

+ {d-2)d^OB^.{x)V with B = {D,C}, (10) 

is decomposed in operators that can be easily classified according to their renormalization prop- 
erties. Here the weight function reads w^ = 1 and w'" = 2x- for scale and special conformal 
transformations, respectively. We introduced the following set of type A and B operators 

Oa{x) = \ {Gi:f , Ob{x) = 5^^^^ (c."9,5;) , Os,{x) = 6^^'^ {oo'^B^) , (11) 

as well as class C equation-of-motion operators 

^o{x) = B;J§,, ^^.(-) = ^V^ + V^g, nM=^'§„ ^d{x)=D^§,. (12) 

The renormalization of Eq. (^) is straightforward and the renormalization of the operator prod- 
ucts and the resulting renormalized conformal Ward identities are given in Ref. |^. 

As a side remark let us note that in spite of the fact that the conformal field transformation 
laws for the dimensionally reduced, from d = A — 2e (and e < 0) to 4 dimensions, theory differs 
from the ones in dimensional regularization by the presence of e-scalar contributions, e.g.0 S^B"^ = 
^v^liu ~ Bp. ^^^ ^~B^ = {2x_xc — x'^ri'DG'^^^ — 2x^B°^ for the gauge covariant variations of four- 
dimensional fields. Nevertheless, the final result for the variation of the action takes the same 
form as in Eq. (|10D but with boson fields being four- dimensional instead. 

3 Superconformal anomalies. 

In four-dimensional space-time the classical action of the A/" = 1 SU{Nc) super- Yang-Mills theory 
in the Wess-Zumino gauge, i.e. 

5,1 ^Jd'xC^ix) = jd'x {-i (g;,)' + '-rr'^' + \ [D^f] , (13) 

contains the Yang-Mills field strength G^^, = d^.B'^ — d^B^ + gf^^^B^^B^^ the Majorana field if)"" 
satisfying the conventional condition ip'^C^'^^ = ip, and an auxiliary field D"^. It is invariant 



■^This choice is legitimate since the infinitesimal conformal variation is linear in d^ and thus does not affect 
the Ward identities, since the anomalous part will show up as a renormalization countcrterm of the product of 
conformal and equation-of-motion operators. 

''The indices fi, jl and /2 refer to the 4, d and 2e dimensional spaces, respectively . 



under transformation of the superconformal group which consists besides the conformal group 
also of the translational, Q, and conformal, S, supertransformations. The latter two are defined 
infinitesimally by their action on field operators as 

S^r = \G%<y,uC - ^/^^sC, 5^Bl = -iCl.r. S^D^ = C T%^\ (14) 

where ( = (^q — i^C,i is the Grassmann parameter. For (^i = (F = Q) we have restricted super- 
transformations, while for C^q = {] [F = S) these equations define the superconformal variations. 

The superconformal group is defined by its algebra from which we will be interested in one 
particular commutator 

[Q,/C_]_=7-5, (15) 

with Q [S) being super (conformal) generators. Note, however, that for the short supermultiplet 
(S" '?/'", D") this commutation relation is modified for action on fermions. To restore it one has 
to use Jackiw's gauge covariant conformal transformation 5^ = 5^ + (jsauge jg]^ where the gauge 
transformation is defined with field-dependent parameter ej^ = (2x^Xy — x'^gfj,i^)B^, instead of the 
conventional 5^ variation defined above. For the action on a space spanned by gauge invariant 
operators this modification is irrelevant. The commutator (p!5|), when applied on a Green function 
with conformal operator insertion, will provide in the supersymmetric limit (identifying colour 
factors) non-trivial relations between the afore mentioned QCD special conformal anomalies. 



To quantize the theory described by the action (|T3D we have to add a gauge fixing and a ghost 
term. We do it via the covariant gauge fixing 

^gf = Jd'xl^-y^{d,B;r + d.co'^V^'u;'^ . (16) 

Although it explicitly breaks the supersymmetry on Lagrangian level, it will not affect gauge 



invariant quantities since the supersymmetry variations ([T^) commute with BRST transformations 
on the gauge fixing function^, i.e. [5^, (5^^^^]_((9^i?^) = 0. 

Translational and conformal supervariation of the action S = Sd + Sgi regularized by means 
of the dimensional regularization, 4 —>■ d = A — 26, leads to 

^S'^S = - {(oOs^) - Og^s^ (17) 

^5'S = ^ {C,A) + (003,) - Or\ (18) 



^We write for brevity (J^i^ST j^g^ga^j of (5^^^'^/(5A, i.e. after transformation the infinitesimal Grassmann variable 
is canceled from the right. We recall that the BRST transformations are given by the set of equations (jBRST^a _ 



where O = J dfixO{x) and the operator insertions read 

A{x) = a,,Gl,r, (19) 

c?3^(^) = 4/'^"' (^"7mV^') (7mV^'^) , c)3'^(^) = ^r'' {ri.i^') {^i,r) , (21) 

with the operator A being the superconformal anomaly [0, ^ in the trace of the supersymmetric 
current, i.e. Qp = ^G'lj^^a^u'-ypip"'. We used in Eq. (^) the identity 

^F^BRST ^^-d,B;) = 2,5^^ST^F (^^-d^B;) + 6^Q^. 

Note that the three-fermion operators O^^ and O^^ vanish in four dimensions by means of Fierz 
rearrangement. Moreover, O^^ can be generated by a special conformal variation of the operator 
©3^, namely, 5^0^^ = i'y^O^^ + 0{e^), where at one-loop level we can safely neglect the remainder. 
Later on we will concentrate on the use of the ^ = —3 gauge in the derivation of the constraints, 
which ensures renormalized supersymmetry at one loop order [Q. This means that the quark 
and gluon anomalous dimensions are equal 7^ = 7^ = 7^. Furthermore, at one-loop order it was 
found that anomaly A does not acquire gauge variant counterterms |]1T|, provided one uses this 



particular value of the gauge fixing parameter. Therefore, we can write to this accuracy eA in 
terms of the renormalized operator 

^A=^[A] + Oia',), (22) 

with d-dimensional /3-function f]^ = —eg + /?. Finally, we write the superconformal variation of 
the action to one-loop accuracy as 

zS'S = ^0[^] + C^O^^ - Of s^ + 0{al). (23) 

Let us point out a further consequence of the ^ = — 3 gauge, which also leads to the equality 
of anomalous dimensions 7^ and /^-function, namely P/g = 7^ = — Q;s/(47r)3A^c- Consequently, 
the renormalization of the conformal variation of the action will be simplified and its integrand in 
Eq. (^) reads in one-loop approximation (see e.g. Ref. [0) 

^ {Oa{x) + Ob{x) + f2o{x) - n^^x) - f2D{x) 

2 

+ 2j^n^{x) + (d - 2) d^PB^^ix)], (24) 



-^ ([OAix)] - n^4x)) - ^^([OBix)] + na^ix) - noix)) 



with ghost anomalous dimension 7^. 



4 Superconformal Ward identities. 

In order to derive Ward identities we have to know the change of conformal operators under the 
superconformal symmetry. Using the rules in (|l^) one finds that the translational supersymmetry 
transformation laws are given by (here and everywhere aj = ^[1 — (~1)"']) 

5« Sji = a, CoV.-ii, (5« si = a, CqV,,, 5« Vji = a,+i CoW,_iz, (5« V} = a,+i CoW,7, (25) 
5^ Vj-ii-i = -7-Co {Sji + S]_,i} - 7-75C0 {Vji + V]_,i} . (26) 

These equations make our comment about the formation of the conformal operators into the chiral 
superfield apparent: Eqs. ( P^IPSJ ) are in one-to-one correspondence with the super symmetric rules 
for the Wess-Zumino multiplet ^. Under superconformal variations conformal operators behave 
as follows 

5'Sl = -a,{j + l + 3) Ci7+V,-n-i, 5'Sl = -a,{l - j) Ci7+V,7-i, 

S'K = -^j+iU + ^ + 3) Cil+Kj-ii-i, 5^Vl = -a,+i(/ - j) Ci7+W,,_i, (27) 

5' V,.u-i = -7-7+Ci {{I - j)S]i^i + {J + 1 + '^)SUi-i] 

-7^757+Ci {{I - 3)V]i^i + {3 + 1 + 2)^|-i^-i} • (28) 



Note, that Eqs. (p5| , p7D do not require Fierz rearrangement and, therefore, they do not change 
their form when the theory is regularized via dimensional regularization. Finally, let us recall that 
the transformation laws of conformal operators under scaling and special conformal variations are 
given by 

5D^O,i = -{l + d{f])) ""O.u ^"i'^O.i = ^a^l{n)''0,l^^. (29) 

Here d{B) = d{G) = d{Q) = 3 and d{F) = | as well as aji{B) = aji{Q) = aji(G) = a{j, 1, 1, 1) and 
aji{F) = a{j, I, 2, 1) with a(j, /, Ui, 1/3) = 2(j -l){j + l + Ui + iy2 + l), where u^ = d^ + s^-1 and s^ 
being the spin of the field 0. Again the scale dimension d^ is chosen to coincide to its canonical 
value in four dimensions. 

Due to difficulties to preserve the supersymmetry of the theory with quantization and regular- 
ization procedures our modest goal will be, therefore, a derivation of the constraints for the special 
conformal anomalies of the QCD conformal operators stemming from the commutator equation 
(|15|) at one loop level only. We will choose the covariant gauge with ^ = — 3, which gives us the 
advantages mentioned above. 

The dilatation and special conformal Ward identities for a conformal operator ^O, which is 
either bosonic {O = B) or fermionic {f2 = F) one, look now very simple (cf. Ref. 0) 

{[''0,i]6^X) = t{il + d{n)) 1 + 7(^)},, (fO,,]X) + ^UfO.i {Oa - ^^^)]A'), (30) 

A:=0 6 

8 



k=0 & 

-2{^fO,lA:JX), (31) 

with a two-by-two unit matrix 1 = l[2]x[2]- Here 7 and 7^^ are the scale and special conformal 
anomalies. It is well known, the scale Ward identity coincides with the Callan-Symanzik equation. 
Thus, 7 is the conventional anomalous dimension matrix and the combination {I + d)l + ^ is the 
scale dimension matrix of conformal operators. Obviously Oa{x) — Q^^{x) = —2Cci{x) is the 
classical Lagrangian (p!3| ) without auxiliary fields. In Eq. ( pO| , |3TD we have also introduced a new 
conventions for the operator insertions weighed with different functions, i.e. O = J dfix 0{x), 0~ = 
J d'^x 2x_ 0{x) (analogously for equation-of-motion operators) and A^^ = J d'^x 2x„ d^OB^i{x). 
The precise definition of the renormalized operator products is given in the next section. 

Let us turn to the renormalized supersymmetric Ward identities. For definiteness let us con- 
sider parity even S operators. From unrenormalized Ward identity (^, the super conformal varia- 
tions of the action ([T7| , ^3D and the operators (2^,27) we can immediately derive the renormalized 



Q and S supersymmetric Ward identities in one loop approximation: 

fe=0 S 

-{{CiOs^) [^!i]X) + {[S^.P^'^'^X). (33) 

Here the super conformal anomaly rk is the residue of the renormalization constant 

r,k = rfi + Kfi + . . . , (34) 

arisen from the renormalization of the operator product 

(Ci[^]) [s;,] = [{CiA) s;,] + a,Ci7+ E v^jv^.-i] (35) 

fc=0 

and induced by the e term of jS^. Note that since the three-fermion operators vanish in four 
dimensions, their product with the bosonic operators will give a finite contribution at one loop 
order. Similar equations hold for Vji with the replacement of index V hj A and Cj by (Tj+i. 

We have neglected infinite terms in the above Ward identities, since they have to cancel each 
other. It is instructive to discuss this issue in more detail for the supersymmetric Ward identity 
of the two- vector S = (^2) 

{[S,i]6QX) = -{6^[S,i]X) - {z[S,i]{6QS)X). (36) 



The variation of "good" component of the fermion field, ip^ = |7_7+V', entering in X may cause a 
divergency on the l.h.s. since it contains a composite field strength G^v Fortunately, in the light- 
cone gauge the latter can be expressed in terms of elementary vector potential B^ and, therefore, 
the l.h.s. is finite by definition. This gauge, together with the use of dimensional reduction which 
implies 5^S = 0, leads to 

{6^[Sji]X)=&mte. (37) 

Since the renormalization of the composite operators is both gauge and scheme independent at 
leading order, Eq. (|37D holds also true for our choice of scheme. Furthermore, in Eq. (|3^) the 
renormalization of the product [SjijO^^ is finite at leading order because O^^ ~ ^'(e) and cancels 
a pole in one-loop diagrams. Consequently, the product [Sji]Oq^^^ = finite, or if divergent it 
cancels singularities in [Sji]6^iIj mentioned above. 

Now we discuss the consequences of Eq. (^Tl). The components of the renormalized operators 
[Sjk] are defined in terms of unrenormalized ones, constructed though from the renormalized fields 
0, by 

[S^i] = E tr'ZshkSl SI, = Z^'sf\ (38) 

6=1 A:=0 

where S-, is expressed in terms of the bare fields 0*^°^ = J Z^cf) and coupling g'^^^ = ^'^g/ J Z^. 
The anomalous dimension matrix of the vector [S] is defined as usual 



In our scheme the anomalous dimensions 

are expressed by the residues of the Laurent expansion of the Z factors Z = 1 + Z^^^e + 0{£~'^). 
The renormalized fermionic operator is defined by the same equation (PB]), however, with the 
2 X 2-matrix Zs being replaced by the numbers Zy. 
From Eq. (^) we conclude that 

k=Ok'=0 \ 1 ■^SJjk 1 ^Sfjk ) \ l^V ]kk' ) 

where we implied that Z^k = for k > j. Substituting the Laurent series into this result, the 
--poles have to cancel. This is ensured by the relations 

t C'^™},. ^^ t^^-i''] + t {''Z^},, ^^ \M = ^. E {4'}, , , [Vki], (42) 

fc=0 ■' k=0 ■' k=0 ■' ' 

t {"41, -. [Vk^.A + 1 {-41, -. [M = -. E {41, [Vh], 

fc=0 fc=0 A:=0 

10 



where [Vki] are independent operators. These when combined together with the analogous results 
for parity odd operators (replace S ^ V, V -^ U, and aj -^ ^j-i) and Zy = Z^ provide 
constraints for anomalous dimensions (when differentiated w.r.t. Ing) |^, |^: 



,.2m 



I '^ J 2n+1.2m+l L ' J 2n,2m L '^ J 2n.2 

{^^4H ={''zm ={zm 

I ' J 2n,2m L "^ J 2n-1.2m-l L '^ J 2n-l,2m-l ' 

{''zm =MH =141} 

L "^ J 2n+l,2m+l L ' J 2n,2m+2 L "^ J 2n,2m+l ' 

I^^ZW} =r^4^| =14' I • (43) 

L "^ J 2n+l,2m+l L ' J 2n+2,2m L "^ J 2n+l,2m ^ ' 

Now let us turn to the renormalization of the superconformal Ward identities. An important 
consequence of these constraints is that the operators are multiplicatively renormaliziable in the 
one-loop approximation, e.g. {^^^l | = |2i^[i] } = q and {^^Z™ } .^^ .^^ = {^^zg' } = {z™ } . 
Thus, in this approximation the classical transformation laws (|27|) for superconformal variation 
remain true also for the renormalized operators. Consequently, the superconformal variation of 
the renormalized operator provide finite Green functions 

{6'[S,i]X) = {[S^S,i]X) = finite. (44) 

Since Cj^ vanishes in four dimensional space-time, ((CiC^3v>) ['^iJ'^) i^ finite. Thus, also the 
superconformal Ward identity (|33D is renormalized up to possible divergencies on the l.h.s. that 
are canceled by the renormalization with BRST-exact operators on the r.h.s. of this Ward identity. 
Note that in our leading order approximation the anomalous term proportional to the f3 function 
is given by a tree approximation. 

5 Constraint equalities for conformal anomalies. 

Having derived Ward identities we are now able to discuss consequences of the superconformal 
algebra. To demonstrate the method, we derive at first the set of relations (^31) between the 
anomalous dimensions arising from the commutator of super and scaling variations 

[<5Q,<5^]_ = 1<5Q, (45) 

which is deduced from the commutator algebra [Q, !?]_ = |Q. Next we deal in the same conceptual 
manner with the commutator (|l^) of super and special conformal transformations, written in a 
symbolical form as 

[5«,5^]_ = -z7_5^. (46) 
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This provides us the desired constraints for the special conformal anomahes of the conformal 
operators. We mostly concentrate on the even parity sector and just state the results for the odd 
one. 

Before we start, let us argue that the three-fermion operators and BRST-exact operators 
will not contribute to the constraints. As already mentioned, the product of the three-fermion 
operators with composite operators provide in one- loop approximation a finite part, which could 
possibly contribute to the constraints. However, its evaluation gives a result that depends on the 
75 and tr 1 prescriptions. Let us explain this point in more detail. Calculating the contributions 
of this operator product one deals with a quark loop that due to different Wick pairings has three 
terms: two of them contain a string of Dirac matrices while the last one is a trace. Obviously, 
in four dimensions the sum gives zero result, recall that O'^'^ and 0_ vanish by means of Fierz 
rearrangement. However, because of e~^ pole in the loop momentum integration, we have to keep 
^-contributions from the spinor algebra that cancel this pole. Obviously, the e part is ambiguous 
for axial channel: the result depends on the handling of 75 in the string of Dirac matrices as well 
as in the trace. Next, since one of the three contributions is given by a trace of Dirac matrices, the 
result depends on the prescription for the trace of the unit matrix tr 1. One of the choices made 
in most calculations is to adopt a fiction for rf- dimensional gamma matrices that still tr 1 = 4. 
However, in those computations the trace appears as a single overall factor and the above choice 
is permissible. It results into scheme dependence for finite part of e.g. one-loop diagrams. In our 
case, since we have an additive trace contribution, we have to continue the Clifford algebra in 
d dimensional space as well, this results into the rule tr 1 = 2['^/^l. This convention produces a 
term involving In 2 reflecting scheme dependence. On the other hand we certainly know that in 
leading order the conformal anomalies do not depend on these ambiguities. So we conclude that 
the contribution of the operator products in questions can not affect the constraints. 

Now we come to the operator products containing unphysical BRST operators. Of course, 
one expects that these operators do not contribute to the physical sector, however, they may be 
responsible for the cancellation of unphysical pieces appearing in the renormalization of products 
containing only gauge invariant operators. On the very end we are interested on relations for 
physical quantities and as we already know from the constraints on anomalous dimensions these 
operator products have to be die out in Eq. (^). For the commutator (|4^), we have the super- 
conformal anomaly and the only difficulty could be a gauge dependent term that is cancelled by 
the operator products in questions. From our previous experience in Ref. |]^, 0, we expect that 
such a contribution is absent and this will be shown by explicit calculation in section ^. So it is 
justified to neglect the whole unphysical sector from the very beginning. 
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5.1 Cominutator constraints for anomalous dimensions. 



First let us demonstrate the derivation of the relations (^) for anomalous dimensions at leading 
order from the commutator of scale and supersymmetric variations, given in Eq. (^Sj), applied to 
the Green functions of composite operators 

(47) 

) times ^. To 



The r.h.s. of this equality is obviously given by supersymmetric Ward identity 
calculate the l.h.s. of the commutator we employ the Ward identities ( PUD and ( 
exact operators being omitted, and find the following contributions 



with BRST- 



([«5. 



ji 



^X) = -^(z [(5«5,,) {Oa - f2^^)] X) + {l + I) ((CoOsv,) [S,,]X) 



g 



-XoEiO^Dl^' ,, 



fc=0 




{{Si,]S''S'^X) 



I, 

g 



SjiiO 



a 



'tp'(p 



6jk 



fc=o I yu Ojk 






(4J 



X). (49) 



To derive the r.h.s. of Eq. 



with 5^0^^ = —^O^^ and the variation 



we have used for ((CoC's^) [Sji\5^X) the scaling Ward identities 
]) of conformal operators. The variation of the action 
proportional to e is neglected, since this Green function is already finite at one-loop order. As it 
was already expected from our previous result ( ^3] ) and a discussion about scheme dependence, 
the three-fermion operator insertion will not affect the constraints at one loop level. 

Now we come to the terms in the above equations proportional to the (3 function. To derive 
we have used an equation (modulo infinite constants which again do not affect 



(50) 



the r.h.s. of (g8 

the constraints, since the latter are basically relations between finite contributions) 

{t[S,i {Oa - ^^^)]X) = {z[S,i] {[Oa] - ^^^) X) - 2{z[S,i]X), 



to get rid of the variation sign on the field monomial 6^X. It results from the study of the 
differential vertex operator insertions in the Green function {[Oji]X) with bosonic conformal 
operator O = ( g^ ) ? and one finds [|^ that the renormalization constant of the operator product 

[O,i][0a], 



3 



ipA{x)][0,,] = i[0a{x)O,,] - 6^'\x) E {Za},, [O,,] - -d^S^'\x) ^ {^a},, [O 

k=0 k=0 



kl-W 



dg 



2e^l Blix) 



di 



SB-{x) 



(51) 
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contains a finite contribution (second term on the r.h.s.) 

The constant X is related to the charge and gluon wave function renormahzation constants by 
the relation X = Zg^fZ^. For the \Sj]\Q^^ we have to use the identity 

i\o^^n^{x) = t[o,in^{x)] - ^(x)j^[o,i]. (53) 

where it is obvious that 

{[0,if2^{x)]X) = ^{[0,l]<|>{x)J^^X). (54) 

So that in Eq. (|50|) the finite piece appears as a consequence of two contributions: ^O part from 
the [Sji] [Oa] product due to finite part in Eqs. (|51|j5^ ) and '^O part from the [Sji] Q^^ by means 



of Eq. (pS]). Similarly, we have for the fermion operators a finite contribution 

i[V,i][Oa] = i[V,iOa] + [V,,] + O(e-0, (55) 

for the trace anomaly, and 

i[V,i]Q^^ = i[V,iQ^^] - [V,i] + 0{e~^), (56) 

for the equation-of-motion insertion. In both cases one should note the factors of one in front of 
the second term on the r.h.s. not 2 as for the quark and gluon operators. Although Q^ = Q^ for 
Majorana fermions (recall ip^'Ju'ip"' = due to the Majorana fiip properties), tp and ip are treated 
as independent variables in the functional integral. Finally, we have 

(25«[5,, {Oa - ^^^)]X) = (z[(5«5,,) {Oa - ^^^)]X), (57) 

which is almost a trivial result. To derive it one uses (pOD and observe that the variation of the 
last term {[Sji]X) in it cancels with the second terms in Eqs. (|55|j56|) so that we are left with the 



r.h.s. of (0). All other terms in the commutator of Ward identities are relatively straightforward 
to handle. Obviously, with the equality (|^) the contributions proportional to the /? function in 
Eqs. (|48D and (^) cancel each other in the commutator relation. 

Subtracting Eq. (^) from (^) and comparing the difference with the supersymmetric Ward 
identity (^) we get, after extraction of independent combinations and identifying both parity 
sectors for fermionic operators, the known supersymmetric relations [^ (see also []2[): 

11 -5 _ 22 'P _ V ^ 

72n+l,2m+l ~ 72ra,2m ~ 72n,2m5 TTl ^ U, 

12 S _ 21 'P _ V <r _ 1 

72n.+l,2m+l ~ 72n,2m+2 ~ l2n,2m+li TTl ^ Tl i, 

21 "5 _ 12 ''^ _ V ^ (r\S.\ 
72n+l,2m+l ~ 72n+2,2m ~ 72n+l,2m5 171 ^ Tl , [00) 

22 S _ 11 ^ _ V ^ 

72n+l,2m+l ~ 72n+2,2m+2 ~ 72n+l,2m+l5 m ^ U , 

/2n+l,2n+l ~ "-' 5 /2n,2n ~ "-'• 
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Obviously, these are the same relations as given in Eq. (^). If one relies on a supersymmetry 
preserving scheme, these equations can be derived in any order of perturbation theory and they 
have been checked at two- loop orderQ 0. 

5.2 Commutator constraints for special conformal anomalies. 

The constraints for the special conformal anomalies of the conformal operators result from the 
commutator (|46|) of supersymmetric and special conformal variations applied to the Green func- 
tion: 

{[S,i]([6^,6^]x)) = -z^^{[S,,]6'X). (59) 



The derivation runs along the same line as above up to the appearance of the superconformal 
anomaly on the r.h.s., given by the Ward identity (^3|), and the absence of finite contributions to 
the renormalization of the product [Oj7][(9^]. Again omitting the BRST-exact operator insertions, 
the commutator of the l.h.s. is given by the two equations 



g 






g 



s n \ / 11 c,<> 12 c,o \ 1 / r-v 5 



^Co E ^M «.Ki?) ' + % ''\\\{\'^-''-''\X). (61) 



The conformal variation of the Green function with three-fermion operator appearing in Eq. (|& 
can be calculated at tree level: 

([(CoOs^) S,i]5^X) = -m,,(B)([(Co03^) S,,i^,]X) - ([(Co5-03^) S,i]X). 

Here we again neglected the BRST-exact operator that arises from the conformal variation of 
the action. Taking into account the Ward identity (33) and the relation 6^0^^ = i'y^O^^, we 
observe again that the three-fermion operator contributions cancel each other in the commutator 
constraint. 



^Here we evaluated the rotation matrices from the conventional MS, in which all next-to-leading anomalous 
dimensions are available, to dimensional reduction scheme. As mentioned in the introduction this procedure does 
not completely fix the non-diagonal part of the rotation matrices. 
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Let us now consider the operator product proportional to the 13 function, which is defined by 



m 



o,i[o 



i''4>) 



6 



6 



X) = {t[OMO^ - (2rJ]X)-{{ d''x2x. U— + ^-^ [0,i] X). (62) 



The supersymmetric variation of this product gives 

(z5« [5,, {OJ - Q-^^)] X) = {^[{S<^S,^) {Oj - f2^;)]X) + z([(Co7-^) S,,]X) 



+ {(j d'x2x.^j^ 5«[5,,]') X) - {6^ (jd^x2x^ L±. + ^l^ is,i]\ X). (63) 

Here we have used 5^ (O^ — ^j^^j = (^07--^+ • • •; with elhpsis standing for the Ca^ operator which 
again is irrelevant in one-loop approximation, since the whole contribution (^) is multiplied by 
the /3-function and thus starts from a.^- The first line of this equation ensures the cancellation with 
the same terms appearing in Eqs. (^,^), however, the remaining equation of motion operators 
in the second line will contribute to the constraints for the special conformal anomaly. We can 
absorb these additional pieces by a redefinition of the special conformal anomaly matrix. The 
action of the equation-of-motion operators is 



d'^x 2x. 



I 

r- j 






, _.^ „,..,. ,,.,,_„ (64) 

^^ fc=0 

where bjk{B) matrix can be found in P, ^ and the fermionic one is evaluated in Appendix ^ Now 
we shift the QQ-entry by ^^jf. -^ ^^Bj'f^ = '^'^ji^ + 2-bjk{B) and set for the remaining channels 
"^^Pji^ = ^^jjk- This redefinition treats the P term equivalently for quarks and gluons and instead 
of Eq. (^ we have at leading order 

The new special conformal anomaly matrix 



ijk- 



(65) 



llT-iC 127-iC 

217-ic 227-ic 



(66) 



for the Sji operators can be found from the conventional quark and gluon ones by the transfor- 
mation 



k jfc 

1 127-10 

fc+1 ^ jk 
1 21pc 
k ^jk 

1 227-'c / 
fc+1 jk / 



I 



1 



2A; + 3 



1 fc+a 



-1 

J+3 

2±3 
i+1 



fc+3 
6 

k 
6 

(fc+3)(j+3) 
6(i+l) 
fc(J+3) 
6(i+l) 



6 

i 

_6 

i 

6 

i+1 



fc+3 \ / QQpc \ 

jk 

QGpc 
J- jk 

GQpc 
J-jk 



k 

J 

fc+3 

i+1 

6 fc , 

i+1 j+1 / 



(67) 
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An analogous convention is introduced for the special conformal anomaly matrix of the fermionic 
operator^ 

r^,{F)^j^,{F) + 2%{F). (68) 



If one compares now both sides of Eq. ( [59| ) in terms of these new conventions we find that 
besides the special conformal anomaly r"^ only the superconformal anomaly, 



'^A],^-2\f, (69) 



jk — jk 1 

arising from the renormalization of [v4.]5j7 in Eq. (|35D, contributes to the desired constraints: 

k=0 fc=0 

-^j E {^Ukin - K} [^^''-i] = 0' (70) 

fe=0 
k=0 k=0 

-^. E {m^) - "^Jk} i^ki-i] = 0, (71) 

fc=0 

and for Vji they read 



fc=0 fc=0 



-^.+1 E {^j-iA^) - K} i^^i-i] = 0' (72) 



fc=0 

E ^^r,^^(5)^.+i[WH-i] + E ''^^(i?)^;^+i[w.-i,z-i] 

fc=0 fc=0 

-^.+1 E {^^(^) - '^/4 t^^'-i] = 0' (73) 

fc=0 

where we have implied r'^^ = for k > j. Extracting the independent components from Eqs. 
( |70| - [73D we finally obtain four equalities for the non-diagonal elements 



-'2n+l,2m+l(-D) + ^2n+l,2m = -'2n,2m(-D) + ^27i,2m = ^2n,2m{^)^ ('4) 

^ 2n+l,2m+l{B) + ^2n+l,2m+l ~ ^ 2n+2,2m+2\B) + ^2n+2,2m+l ~ ^ 2n+l ,2m+l\^ ) i l'^) 

^ 2n+l,2m+l\B) + ^2n+l,2m+l = ^ 2n,2m+2\B) + ^2n,2m+l = -^ 2n.,2m+l (-^ )) ('") 

-f 2n+l,2m+l(-^) + ^2n+l,2m = ^ 2n+2,2m,\B) + ^2n+2,2m — ^ 2n+l,2m\^ )i ('') 



^ Here the fermionic conformal anomaly "i'^{F) has the same structure as the quark one, namely, 7'^(i^) = 
—b{F)j{F) + w{F), with anomalous dimension 7(F) of the quark-gluon operator V (U) and a part w(F) deduced 
from the renormalization of V (Z^) with the trace of the energy-momentum tensor. 
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Figure 1: One-loop diagrams which give rise to divergences in the product of the renormahzed 
operator insertions ^[^jpOj;]. 

with n > m, and six equations for the diagonal elements, which are of no relevance in prediction 
(0) for the anomalous dimension matrix. Eqs. ( [71[ - [77| ) are the main results of this section. The rest 
of the paper is devoted to a consistency check of our results for the special conformal anomalies 
by evaluating the Z\-anomalies at one-loop order. 



6 Evaluation of superconformal anomaly. 

To simplify the calculation of the renormalization constants A, it is advantageous to use the 
light-cone position formalism and to this end we introduce the non-local light-ray operators 






1 T„ 



r4^^2n)r$'''[K2,Ki]^UK,n), 






(75 



The calculations of the mixing kernels for the product [A] [O] with fermionic operators are straight- 
forward and go along the line of Ref. |13[ (see also appendix B of Ref. |^ for a recent review of this 
formalism in the case when the non- forwardness is essential). The diagrams are represented in 
Figs. and ^ for quark and gluon operators, respectively. Since these renormalization constants 
must be gauge-independent, we have performed the computation with a general ^ and indeed 
observed its cancelation in the final result. This shows that the BRST-exact operators can not 
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Figure 2: Same as in Fig. ^but for the operator product i[^]pOj7]. 
contribute to the constraints. The generic form of the result is 



X (^Ci7+ O'j {xKi + XK2, yni + y^a) =F (fi^i ^ k,2), (79) 

where the — (+)-sign in the second term corresponds to the parity even quarks and parity odd 
gluons (parity even gluons and parity odd quarks) and we use throughout the convention x = 1—x. 
Skipping details (see Appendix |B]) we give the result 






S{x) 



Siy) - 2 + S{y) + (1 - y)S{x) 



«7^^(x,y) = «7^^(x,l/), (80) 

^7^^(a;,y)=^7^^(x,y)+%, (81) 



with the +-prescription defined conventionally by 

1 



X 



i-5(x)r^. 

X JO y 



We have kept in the quark sector the gluon momentum ^2+ unintegrated, which stems from the 
exponential in the Feynman rules for the operator insertion Vertex x exp {—inikij^ — iK,2k2+)- We 
merely substitute it by the corresponding parton momentum fraction when passing to exclusive 
type kernels (see Appendix 0). Note that the result for the polarized quark kernel is the same 
as for the parity even case. This could be anticipated, since for divergent parts we can use 
anticommutativity of 75. However, there is a sign change in the contribution of the diagram in 
Fig. ID (c) but it gives the contribution ±x6{l — x — y) with +(— ) sign for even (odd) parity 



19 



and together with x6{l — x — y) from the Fig. |1| (d) gives a vanishing result when the symmetry 
property of the corresponding quark string operators is used. For gluons the difference between 
the parity even and odd sectors comes entirely from the contributions of the graphs in Figs. § (b), 
(c) and (e). The contributions from (d) and (g) annihilate each other, while remaining diagrams 
give identical^ results for V and A channels. 

Evaluation of the conformal moments according to the method spelled in Appendix gives 
for gluons 

|iV.)" ^7^W = ^%-i,.| ^^ + 'f^+'^^' + {k + 2){k + m,, + {-ir\k + 2)ik + 3) -.}, 

(82) 
with — (+) sign for vector (axial) channel; and for quarks 

^iV,)"'«7^W = 2^|^|(-l)^+'=(A:+2)-l-(-l)^+'=(j + l)(j+2)^,,-(j + l)(j+2)r,,|. (83) 

To simplify the presentation we have introduced the matrices, depending only on logarithmic 
derivative of the Euler integral ipi^) = ^l^-^(^)? ^^^ 

^,k = i^(^^^^]-i^ij + k + A)-i;(^)+i;ij-k), 

^jk = ^{j + k + 4)+^{j-k)-ij{k + 4)-^{k + l), 

Sjk = ij{j + k + 4)+ij{j-k)-2ij {j + 2) . (84) 

The conformal moments TZjk of the kernels TZ{x, y) are related to those in Eq. (^) by a normal- 
ization factor and, therefore, the super conformal anomaly is 

"Aik = f^cu^ ^7^,•fe + ^cj" ^7^,•fe) with a = 1, 2. (85) 

We should note that for even and odd j — k the ■j/'-functions with argument depending on both j 
and k enter only in the particular combination 

Ajk = ^ {^-^^\ - ^ (^) + mj -k)- Hj + 2) - ^(1), 

which arises in all special conformal anomalies 7*^, too |^. 

The final step is to insert our findings for the superconformal anomalies given in Eqs. ( p2| - p5| ) 
and our results for the eight special conformal anomaly matrices from Ref. 0, rotated to '*^r''^-basis 
by means of (|67|) , into the four conformal constraints (|7^[77| ). Indeed, we find that all of them are 
identically fulfilled. Thus, the universality of the special conformal anomaly matrix for even and 
odd parity sectors arises from the form of superconformal anomaly in A/" = 1 Yang-Mills theory. 

^Discarding the fact that contributions with ki <-> K2 enter with opposite signs. 
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7 Conclusions. 

The use of A/" = 1 supersymmetry allows to find non-trivial relations for both scale and special 
conformal anomalies of QCD conformal operators. In the former case, the constraints, at leading 
order in dimensional regularization scheme as well as to all orders in a supersymmetry preserving 
scheme, involve only anomalous dimensions. There are six constraints for the diagonal and four 
for the non-diagonal elements of the four entries in the anomalous dimension matrices in the 
parity even and odd case. Although, in the supersymmetric limit the colour factors have to be 
identified, the constraints have predictive power also beyond the leading order approximation. 
As it has been observed in Ref. [0| there is a subtlety in finding the supersymmetry preserving 
regularization scheme in the non-forward case. Standard dimensional reduction [§ does not serve 
its purpose. However, an existence of a supersymmetry preserving scheme was proven by the 
ability to fulfill the constraints by a multiplicative renormalization of anomalous dimensions. 

There exist four constraints for the special conformal anomalies. They contain new ingredients, 
Z\, due to the superconformal symmetry breaking by the trace anomaly. A, in the spinor current. 
The four symmetry violating entries arise as a counterterm in the product of A and conformal 
operators. Thus, there is no predictive power in these constraints, however, they serve for a 
consistency check of the special conformal anomalies. 

Let us mention for completeness, the situation with the two further leading twist-two conformal 
operators appearing in the definition of the so-called transversity distributions. It is well known 
that these quark and gluon operators do not mix with each other under renormalization due 
to different spin representations w.r.t. Lorentz group. Their forward anomalous dimensions are 
related in the supersymmetric limit by one constraint. However, unfortunately, there exist no 
constraints for the non-diagonal elements of conformal anomalies that do not involve anomalous 
contributions of fermionic operators ||14|| . 

We have evaluated the superconformal anomalies in one-loop approximation and found that 
our relations are indeed fulfilled with the known results for the special conformal anomalies from 
Refs. [^ 1^ . In other words, we can now reconstruct the special conformal anomalies of the gluon- 
gluon and gluon-quark sector from those in the remaining two channels, provided the former two 
would not be already known. Furthermore, we know that the latter give us the rotation matrix to 
the conformal scheme in which the conformal operator product expansion of two electromagnetic 
currents is valid. Since the Wilson coefficients of this expansion in both of these remaining 
channels are fixed by the ones of the forward kinematics and coincide with the explicit calculation 
in the minimal subtraction scheme |T^, rotated to the conformal scheme fl^, we have a complete 



consistency check of all special conformal anomaly matrices evaluated at one- loop level. Since these 
conformal anomalies induce, together with terms proportional to the (3 function, the non-diagonal 
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elements of the anomalous dimension matrix at two-loop level, we have also a complete but indirect 
check on their correctness from the superconstraints (^). Of course, in the flavour non-singlet 
sector the anomalous dimensions arising from our prediction coincide with the conformal moments 
of evolution kernel calculated at two-loop level |1^ . 



Altogether, we have a complete list of consistency checks for the field theoretical treatment 
of conformal anomalies, their evaluation at one-loop order and next-to-leading predictions arising 



from their use. This also supports our results for the transversity sector [|1^], which has been 
treated in the same manner. 
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A Evaluation of the fermionic 6-inatrix. 



The fermionic 6-matrix is defined by Eq. (|6^ ) . In order to evaluate the matrix it proves convenient 
to make the Fourier transform on the fields, i.e. (f){x) = J d'^k e~'''^(j){k) and introduce a set of new 
variable of the integrand: Y = k2+ + ^1+ and X = ^^"'"~ ^+ , so we get 

j d''x2x.^lj— V,i = 2iQ, j d''k^d%Y'-^[tp\^'^\x)}G^^{k{)-fj^^{k2) 



2tY,^B,k{2,l)Vui-i, (A.l) 

fc=0 ^^^ 



with differential operator £ = / + (1 — X)^ and Bjk being defined by the integral 



B,.(a,/3) = /' ix'!!^^^Pt'\x)CPr\x). (A.2) 



The weight and normalization factors are given by standard equations 

w{X\a,p) = (l-XY(l + Xf, 



„+/3+i r{k + a + l)r{k + (3+l] 



{2k + a + l3 + l)r{k + l)r{k + a + f3 + 1)' 

Using Rodriga's formula for Jacobi polynomials and integrating k (and k + 1) times by parts we 
come to integrals which can be easily evaluated by means of the result 

/' dX,n(X\a,-r)Pr'(X) = (_iy2°.... n« +i + l)r(/? -7 +m, + 1) ,A3) 

7-1 r{j + i)r{fj-'j)r{a + -f + j + 2) 
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So that we finally obtain 



r{j + a + l)r{k + a + P + l) 



r{k + a + l)r{j + a + p + l)' 
Identifying a = 2 and /5 = 1 we get the result 



where 6 



jk 



{1, ifj>A;; 0, ifj<fc}. 



(A.4) 



(A.5) 



B Renormalization of the operator product f^lfC^]. 



For the calculation of Z-matrices one uses in the light-ray formalism the usual momentum space 
Feynman rules and Vertex x exp (— mifci^ — ii^2k2+) for the non-local operator where 'Vertex' 
stands for Dirac or Lorentz tensor. Introducing the Feynman parameters x, y for the propagators 
we reduce Feynman integrals to the form 



J * TZ{x, y) X Vertex, with J 



dx dye{l-x-y)e-''''+^^^-''^'''+'-'^^^-'''^+^y^'+^^-y'^^^\ (B.l) 
/o Jo 

where the exponential corresponds to the Fourier transform of the coordinate dependence of the 

string operator 'after evolution', C((l — x)ki + xk2, V^i + (1 ~ y)i^2))- To calculate the divergent 

part of the operator product [A\ [O], we have in addition to take the Feynman rule for the anomaly 

[A] defined in Eq. (|1|). 

Our calculation has been performed with arbitrary gauge fixing parameter ^, which has can- 
celed in the sum of diagrams. Therefore, we present the results corresponding to the definition 
( |79D for the separate contributions in ,^ = 1 gauge. 

In the quark parity even and odd sectors we get on the diagram-by-diagram basis from Fig. |I] 



1 



-k2+^n^d) 

-k2+^n^e) 

-k2+'^n^f) 



-{l-y)5{x), 

X 

±-6{l-x-y), 



|)^w-i 



S(x) — S(x), 



6{x) + -i3-y)6{x), 
6{x) + ^ (l + I) 5(^) 



(B.2) 
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where + (— ) sign in the contribution of the diagram (c) stands for even (odd) case. In these 
results we dropped the contributions of the type const ■ 6{x)6{y), since they do not enter into the 
physical part of the constraints ( [f^ - |77D , namely for k < j. In the sum of Eqs. ( p.2|) the term 
I ■ (5(1 — x — y) for vector case and |(1 — 2x) ■ 6{1 —x — y) for the axial one cancels with the ki <-^ K2 
contribution ('+' and '— ' sign, respectively) in Eq. ( |79D and we get the result in Eq. (pOl). 
The gluon case is calculated from the graphs in Fig. ^ and reads for vector channel 

' -2" 



/C2+ /<-(a|*) 



k2+S{y) 
3 



- 2 + 2(5(x) , 



2+ f ' 



-A;2+^7^^ 



3 \ 3 i 

(6) - I 2 - o2/ ) ^2+ - -(1 - 3;)A;i+ + -k ((1 - y)k2+ + a;A;i+) (y/c2+ + (1 - x)ki+) 



+ -(1 - a;)5(|/) {(1 - 'ix)ki+ - 2k2+ - m(l - x)A;i+ (xA;i+ + k2+)} 



-K2+ l<'(c) 



yh 



'l + x)ki+, 



-A;2+^7^J^) = {2 - -y) k2+ - -(1 - x)ki+ + -k {{1 - y)k2+ + xki+) {yk2+ + {l-x)ki+), 



— /C2+ /<-(e) 



-(3-a;)A;i+ 



1-2) ^2+ - 2'^ ((1 - y)^-i+ + ^^1+) (2/^2+ + (1 - x)ki^ 



1 - 5(a;) 



(B.3) 



-A;2+^^(/|$) = fc2+(5(y) 
-fc2+''7^;;l^^) = 5(y) {-(1 - x)A;i+ + 6{x)k2+] 

— K2+ /■C(g) — K2+ /C(^-), 

with K = K,2 ~ 1^1- The subscript $ (-?\^^) stands for the contributions originating from the 
expansion of the path ordered exponential (non-abelian part of the field strength tensor). The 
axial case differs from the previous one only in the contributions of diagrams (b), (c) and (e) which 
are 

7 \ , /3 7 



- A;2+^7^^ 



(fe) 



\ ~ 2 V ^^^ ^ (2 ^ 2 V ^^^ ^ 2'^ '^''^ ~ ^''^^^ ^ ^^^^'' ''^^^^ + (1 ~ x)ki+) 



x 



-k2+''nf 



2 
3 5 



+ 



-h+^T^tc) = -yk2+ - {I - x)k^+ 



6{y) { (1 — Sx)ki+ — 2k2+ — iKx{l — x)kl_^_ — m(l — a;)A;i+A;2+ !• 
— xjki^, 

y]h+' 



(B.4) 



5 \ i 

1 - -y j k2+ - -n ((1 - y)k2+ + xki+) {yk2+ + (1 - x)/cn 



'(') V2 2' 
Summing the separate terms we have to use the formula (see Ref. for a general result) 

J * {k^+ [1 - (1 - x)6{y)] + k2+ [1 - (1 - l/)5(x)]} = 0, (B.5) 

to reduce the result to its final form ([81|). Then we use the equation k2+B^{k2) = iG\^{k2), valid 
to leading order in the coupling, to reconstruct the field strength from the potential. 
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C Conformal moments of regularized kernels. 

The conformal moments of the transition kernels derived in the body of the text are defined 
according to 



n],^ j\t^^^j^^Pi'^'\t) J^x 1^ 'JZ{x,y)q{[l-x-y]t-x + y 



with the weight w{2,l\t) = (1 — t)^(l + t) and normalization nfc(2, 1) = 8 



(fc+i) 



(fe+2)(A:+3)- 



(C.l) 
We have 



for quarks Cj = Cj and C*^ = CjLi for gluons. From this equation it is straightforward to 
evaluate the moments of all parts of the kernels (along the line of Ref. |^) except of the ones with 
+-prescription since in this case we obtain the result in terms of derivatives of hypergeometric 
function w.r.t. its indices which is not easy to handle. In these cases we have to modify our modus 
operandi and develop a more efficient machinery which leads to more tractable expressions. It can 
be achieved according to the re-expansion of the integrand making use of the orthogonality for 
Gegenbauer polynomials. To be more specific let us consider the gluon kernel [l/x]^6{y) which 
in the momentum fraction formalism translates into [6{t — t')/{t — t')]^. We use the following 
regularization of singular distributions 

/■I dx , , , /■! dx 
Jo 



[X 



\l-£ 



[X 



;i)]- 



(C.2) 



[1 — xj+ Jo (1 — X) 

Then using the representation of the Gegenbauer polynomial in terms of hypergeometric function 
2-Fi and using Rodriga's formula for Jacobi polynomials we integrate k times by parts to get 



Mf, 



9{t - 1') 



(-1) 



i+ 



jk 



,(fc + 2)(fc + 3) 

6£ 



^ 12T{j)T{k + 2) Jo 



X 



lE-l 



dyy''+\l - y)''+' 



X- 



d' 



iFi 



3 



xy 



(C.3) 
The simplicity 



d{xyY' 

The first term on the r.h.s. originates from the 0(1) contribution in Eq. 
of the consequent analysis depends on the handling of derivatives acting on 2-^1- If we merely 
differentiate it k times as it appears and perform y integration this gives 3F2. Finally, after x 
integration Eq. ( p.3|) will be proportional to the derivative of 4F3 w.r.t. a low index [see later Eq. 
( |C.12| )]. Fortunately, it is possible to avoid this if one notices that k — 1 ^-differentiations and y 
integration leads to a 2-F1 function with shifted indices. After the last integration is done one ends 
up with an expression with e derivative acting on 3F2. In order to achieve this let us re-expand 
the action of k^^ derivative d/d{xy) on 2-F1 w.r.t. the complete set of polynomials C; (2x — 1) 



d 
dx 



iFi 



-J + k,j + k 

k + 2 



X 



00 

1=0 



-1,1 + 2k + 3 
k + 2 



X 



(C.4) 
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with expansion coefficients 



Cj-kl 



[l-(-ir'i%-fc,m(2/ + 2A; + 3) 



r(/ + 2A; + 3) T{j -k + l) 



r(/ + l) r(j+A; + 3) 
easily obtained from the orthogonahty of polynomials. Then we get 



(C.5) 



M% 



12 r(2fc + 5)r(j-fc + i) 1^^^^ ^^^ + i) + 5^ 



£=0 



^ ^ /-/,/ + 2A; + 3,A;+l 

X > Cj-kl 3-^2 

^ 3 k,l3 2y 2k + 5,k + l+6 



1 . 



(C.6) 



Now the way to handle the derivative of 3F2 is rather straightforward. First we use the fundamental 
identity for 3F2 (/ G IN) 



3-^2 



— /, I + a, P 



r(7)r(7-a) 



-/, I + a,e 



3-^2 



r(7 + /)r(7-a-/)' '\l + a-7,/5 + £ 



(C.7) 



where we substitute 7 = a + 2 — p with p —>■ 0. Then the expansion w.r.t. e is easy to construct 



3-^2 



— /, I + a,e 
p-l,(3 + e 



1 = l + e 



'l{l + a) 



(C.8) 



+r(p -1)1 + (-1)' [(/ - 1)(/ + « + 1) + /5] 



T{l + l + a-(3)T{(3y 
r(2 + a-/3)r(/3 + 



+ 0{e' 



and together with the identity 



we find 



3-^2 



— /, I + a, [3 

a + 2,p + € 



r(p-i) 

T{2-p-l) 

r(a + 2) / 
r(/ + a + 2)1 



l + e 



i Tjl-l + p) 

^ r(2-p) 

l{l + a 



P 



{Si,o + 5i, 



(C.9) 



(CIO) 



Using the results we have just derived, we perform in a last step the summation in Eq. 
according to the formula 



d_ 

We 



^ ^ /-/,/ + 2A; + 3,A: + l 

£=0 z=o 

1 



2k + b,k + l + e 



1 



1 



(j-/c)(j + /c + 3) 



(A; + l)(A; + 2) (A; + 1)(A; + 2)(A; + 3) 

3 ~k' 



(-ly+v i ^^'l^ ] + [1 - {-ly^'Mj + k+4)- (-i)^'+v 



A; + 2 [' '' "V 2 
+ [1 + (-1)^'+V(J - ^) - Hk + 4) - 7/'(l) 



(C.ll) 
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As a by-product we verify the following formula for the derivative of 4F3, which is difficult to 
derive by other means 



r(j - k)V{2k + 5) 'de 



-J + A; + 1, J + A; + 4, A; + 2, A; + 1 



4-^3 



e=0 



2A; + 5, A; + 3, A; + 1 + e 



1 



{j-k){j + k + ^) 



V+k 



{k + l){k + 3) k + 1 

j ~k 



^ 



'j + k + A^ 



+ [l-(-iy+'=]^(j + fc + 4) 



-(-ly+V 



+ [1 + (-iy+^]^(j - A;) - ij{k + 4) - ^(1). 



(C.12) 



Slightly different procedure holds for quarks. In this case due to the presence of the momentum 
fraction k2+ = ^{ki + k2)+ we re-expand the integrand, modified by adding a constant, 



M^k^ 



2 9{t' - 1) 
l-t it' -t) 



;-i) 



jA; + 2)(A; + 3)(j + l)(j + 2) 



+ J jk 



/o Jo 



k + 1 



(C.13) 



xy 



iFi 



-J,J+3 



in the following series 



X 



2 



X 



2F1 



-J,J+3 
2 



0)\ = j:djk2F, 

k=0 



-k,k + 3 
2 



X 



(C.14) 



with the expansion coefficients 

djk = —Oj-i,k 



3 + 2k 



ij-k){j + k + 3). 



Consequent integration and expansion in e requires Eq. ( |C.10| ) as well as the following result 



(C.15) 



3-^2 



— /, I + a, P 
1 + a,(3 + e 

r(i + a) 



(C.16) 



^^o + ^r(/) 



'r(i + a - /3 + /)r(/3) 



4)' ^,,1 + Oie' 



T{l + a + l) [ '" '' \T{l + a-l3)T{f3 + l) 

which can be deduced using the same recipe as presented above. Final summation gives us the 
result in Eq. (Bl). 
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